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Experimental studies were conducted on pressure patterns, gas inges- 
tion, and sealing capacity of viscoseals. Oil, water, and liquid sodium 
were used as the sealed fluids. Results indicated that when the groove- 
land pairs connected directly with the pressurized cavity (without an 
intermediate circumferential groove), an end effect, or ineffective seal 
length, was evidenced by a sharp decay of pressure along the land leading 
edge. Results also disclosed that a helically grooved rotor viscoseal has 
increasing gas ingestion rates with increasing Reynolds number when the 
viscoseal liquid interface becomes unstable. However, if the helical 
grooves are in the housing (smooth rotors), gas ingestion could be elimi- 
nated by increasing Reynolds number. Secondary grooves on viscoseal 
land areas improved sealing capacity and reduced power absorption. So- 
dium was sealed at 300° to 625° F with negligible liquid loss. 

From a theoretical analysis, a set of two-dimensional equations (that 
includes both convective inertia and viscous forces) were formulated for 
a quasi-two-dimensional flow field in parallel groove geometries which 
can be used as a mathematical model for the viscoseal. The mathematical 
model consisted of an infinite smooth flat plate moving relative to a fixed 
infinite parallel groove-ridge surface with the clearance filled with an in- 
compressible, homogeneous, Newtonian fluid. A modified Reynolds num- 
ber evolved that is based on a groove-ridge characteristic length parallel to 
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the smooth plate velocity vector. Application of the modified Reynolds 
number to published experimental results reveals that a significant im- 
provement (previously attributed to onset of turbulence) in bearing num- 
ber or sealing coefficient occurred when the modified Reynolds number 
was near unity; therefore, it is argued that the strong influence of con- 
vective inertia, and not turbulence, is the main cause of this improve- 
ment. For this reason a turbulent flow model should include the convec- 
tive inertia effects. The consideration of convective inertia effects re- 
veals that fluid density p, aspect ratio b/h Q , and number of groove - 
ridge pairs N are additional parameters of pressure generation depen- 
dence that are not revealed by creeping flow analyses. The analysis 
presented herein provides a unified approach which was previously 
covered by the creeping flow and the semiempirical turbulent flow anal- 
yses. No restriction was placed on the net flow; therefore, the analysis 
applies equally to hydrodynamic seals, bearings, and pumps. An arbi- 
trary groove-ridge shape is also permissible in solving equations. Con- 
sideration of convective inertia effects suggests the principle of viscoseal 
axial pressure generation and attainment of zero leakage. 

INTRODUCTION 

The viscoseal is a special case of a general class of parallel groove 
geometries which have many potential applications for use as pressure 
generation devices. For example, the parallel groove has been applied 
in spiral groove thrust bearings (ref. 1), spherical and conical spiral 
groove bearings (ref. 1), helical groove pumps (ref. 2), as well as heli- 
cal groove seals (viscoseal, ref. 3). 
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The viscoseal is potentially useful for sealing liquid metals, such as 
mercury, potassium, and sodium, in space electric-power generation 
systems that require 1 to 3 years of unattended operation and near zero 
leakage (refs. 4 and 5). As pointed out in references 5 and 6, the visco- 
seal has inherent reliability and long life because of the absence of solid 
surfaces in rubbing contact. 

To date, analyses of parallel groove geometry have been limited to 
creeping flow solutions where convective inertia effects have been ne- 
glected. Convective inertia is the fluid mass acceleration due to a spa- 
tial velocity change, for example, a change in velocity due to an obstacle 
in the flow path. This is contrasted to an inertia effect caused by a tem- 
poral velocity change, for example, a variation in velocity at a point in a 
bearing due to shaft runout (which is a function of time). The creeping 
flow analysis is restricted to flows where the modified Reynolds number 
(see DISCUSSION) is much less than 1; that is, 

Re = Re L(“) 2 << 1 
or 

Re 0.01 

where 

Re, = 

L v H 

-< 1 
L 

Creeping flow analyses for pumps, screw extruders, seals, and 
thrust and journal bearings have been published by many authors (e. g. , 
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refs. 1 to 26). Appendix C contains an outline and discussion of the pub- 
lished creeping flow solutions. The main criticism of the creeping flow 
solution (generally used in the laminar flow regime) is that its validity is 
restricted to Reynolds numbers below most engineering applications. 

A number of semiempirical and turbulent flow formulations have been 
devised (refs. 4, 18, 27, and 28) for operation in a regime designated as 
turbulent, but they neglect the convective inertia terms which are impor- 
tant physically. The main criticism of these semiempirical analyses and 
turbulent flow model is that the existence of the turbulent regime is as- 
sumed. Turbulence, however, may not be necessary to explain reported 
experimental results and may not even exist. Even if turbulence does 
exist, convective inertia effects are thought to be far more significant. 
Kettleborough (ref. 29) numerically analyzed the slider bearing with in- 
ertia, turbulent, and viscous terms considered. When inertia only was 
considered, the results were in qualitative agreement with published 
(turbulent-attributed) slider bearing experimental results. Kettleborough 
concluded that the turbulence term did not appear to greatly affect the 
operation of the slider bearing. 

Golubiev (ref. 24) applied centrifugal pump similarity laws to the 
helical groove seal (including operation where the rotor and housing are 
both grooved). This case considers an inviscid fluid (viscous forces ne- 
glected) whereby the pressure generation is found from Bernoulli' s equa- 
tion modified by a geometric constant, where AP varies as the square of 
the velocity. The viscous force cannot be neglected for two reasons: (1) 
the viscous force causes fluid motion, and (2) the viscous force contributes to 
maintaining the generated pressure gradient. Unfortunately, the creeping 
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flow optimum geometry has been extended to be valid for operation beyond 
the creeping flow regime. Since the convective change is the main 
driving force, a complete reevaluation of the optimum geometry including 
the groove -ridge shape must be investigated before conclusions on visco- 
seal performance can be made. Experiments show that the optimum 
geometry changes for flow beyond the creeping flow regime (ref. 27). 

Several authors (refs. 20 and 21) reported that the viscoseal had a 
certain ineffective length at the high-pressure end of the seal. This has 
been called the end effect, and it has been suggested (ref. 20) that this 
ineffective rotor length be substracted from the calculated wetted-seal 
length in order to compare theory and experiment. Muijderman (ref. 1) 
derived theoretical pressure patterns for the spiral groove geometry on 
plane surfaces and developed an analytical method of correcting for the 
end effect. However, pressure patterns in the viscoseal and end-effect 
alterations of these patterns have not been determined experimentally. 
Existing experimental data consists of average pressure measurements 
from which the existence of an end effect is deduced. 

Gas ingestion from the low-pressure to the high-pressure end of the 
seal is reported to have occurred when sealing potassium (ref. 4) and 
when sealing water (refs. 5 and 27). Although this gas ingestion charac- 
teristic may present no problem with a vacuum on the low-pressure end, 
it is highly undesirable in other potential applications because of con- 
tamination and/or detrimental effect on sealing capacity. Reference 24 
reports that in attempting to seal a liquid, the viscoseals (combining 
helical grooves in rotor and housing) worked on an air-liquid emulsion, 
and reference 25 reports that viscoseals (helical grooves on shaft) 
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operate with a mixture of oil and gas. These modes of operation could be 
due to gas ingestion. Leakage observations reported in reference 21 on 
"seal breakdown" and in reference 4 on "secondary leakage" could result 
from gas ingestion into a closed cavity. 

The objectives of these studies were: 

(1) Study experimentally the pressure patterns in the viscoseal with 

particular emphasis on the pattern near the seal end (end effect) 

(2) Investigate experimentally gas ingestion and sealing capacity of 

the viscoseal 

(3) Develop a physical flow model for a theoretical analysis of the 

viscoseal 

The data and analysis contained herein is covered in detail by the 
authors' papers (refs. 5, 6, and 30). 

SYMBOLS 

a ridge or land width 

a’ ridge or land width parallel to viscoseal centerline 
a acceleration vector 

b groove width 

b’ groove width parallel to viscoseal centerline 

C specific heat of fluid 

Cp specific heat of fluid at constant pressure 

c clearance between ridge and smooth flat plate, or radial clearance 
unit normal directed toward axis 
e t unit normal tangent to trajectory 

F x body force in x -direction 

Fy body force in y -direction 
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F_ body force in z -direction 

G dimensionless factor, function of helical geometry 

g acceleration of gravity 

h Q step or edge height (or groove depth) 

k thermal conductivity of fluid 

L characteristic length, length of groove-ridge pair in smooth plate 
relative velocity vector direction 
L ? axial length, normal to plate velocity 
N number of helix starts 

n integer 

AP pressure differential 

P* static pressure at reference state 
P c static cavity pressure 

Pr Prandtl number 

Q net volume flow rate 

Q* heat source energy 

volume flow rate normal to plate relative velocity direction 
Q| volume flow rate in plate relative velocity direction 
q fluid velocity in x-z plane 

R radius of rotor 

Re Reynolds number 

Re* modified or reduced Reynolds number 
T* temperature at reference state 
t time 

U smooth flat -plate velocity or rotor surface speed 

u velocity in x -direction 
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V 

characteristic velocity in the y-direction 

V 

fluid velocity vector 

V p,x 

pressure flow velocity in x-direction 

V p,z 

pressure flow velocity in z-direction 

V 

velocity in y-direction 

w 

velocity in z-direction 

X 

coordinate along ridge -groove 

y 

coordinate across film (between plates) 

Z 

coordinate across ridge-groove 

a 

angle between relative velocity direction and parallel groove- 
ridges or helix angle 

K 

film height 

V 

coordinate normal to plate relative velocity vector 

T " T S 

e 

temperature ratio - D 

T M- T S 

A 

bearing number 

XG 

empirical sealing parameter (defined in ref. 5) 

M 

absolute viscosity of fluid 


kinematic viscosity of fluid 

V* 

kinematic viscosity of fluid at reference state 

« 

coordinate in direction of plate relative velocity vector 

p 

density 

p* 

density at reference state 

a 

shape parameter 

& 

dissipation function 

<P 

angle traversed in polar plane 
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V vector differential operator Del or Nabla 
o 

V Laplacian operator 
Subscripts: 

c based on clearance or film thickness 

g crrnnvp 

L based on characteristic length 

M moving surface 

n normal direction 

r ridge or land 

S stationary surface 

t tangential direction 

z across ridge-groove 

e plate relative velocity direction 

t ) normal to plate velocity direction (axial length) 

APPARATUS AND PROCEDURE 
Experimental Study of Pressure Patterns 
Figure 1 shows a schematic drawing of the experimental apparatus 
and viscoseal assembly used to study end effects and pressure patterns. 
The rotor of the viscoseal is attached to the power input shaft which is 
supported by externally pressurized gas journal bearings. Axial thrust 
due to sealed fluid cavity pressure against the rotor is resisted by the 
externally pressurized gas thrust bearing. The power input shaft is 
driven by a variable -speed electric drive and step-up transmission. A 
magnetic pickup monitors the shaft speed. The viscoseal housing is at- 
tached to a support shaft and roller bearing assembly, which permits 
axial adjustment. A pump pressurizes the seal housing cavity and circu- 
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lates a cooling fluid flow from the reservoir, through a heat exchanger, to 
the test cavity and then back to the reservoir. This cooling flow is neces- 
sary to control fluid film temperature in the viscoseal. The steel rotors 

-fi 

have a centerline average surface finish of 20X10" inch and are mounted 
within 0. 0002 inch of the total indicator reading. The viscoseal housing 
is constructed from acrylic plastic and contains thermocouples and pres- 
sure taps as shown in figure 2. Pressure taps, in a line parallel to the 
axis, are located at groove edges and midland and midgroove positions; 
pressure taps are also located near both edges along one groove. Two 
sets of four pressure taps (90° apart) were used to aline the housing with 
respect to the rotor within 0. 0004 inch as determined by calibration. 
Thermocouples were located flush with the bore and groove root in a line 
parallel to the axis and at each midgroove and midland position. Figure 2 

also gives the groove dimensions used in this evaluation. The housing had 

1 ° 

a 2. 001 -inch bore and a 14-^ helix angle. Groove widths of 0. 19 inch and 
land widths of 0. 13 inch were provided by using five helix grooves (five 
starts). These geometric proportions, which are based on the optimum 
relations given by reference 19, included a 0. 015 inch groove depth and a 
rotor radial clearance of 0. 005 inch. 

The kinematic viscosity-temperature relations for the mineral oil used 
were found to be the same before and after use in the test. These values 
can be found in table I. 

Experimental Study of Gas Ingestion and Sealing Capacity 
In studies on sealing water (fig. 3), the water pressure and flow to the 
pressurized cavity (at the viscoseal high-pressure end) were controlled by 
throttling valves on the inlet and outlet water lines. This arrangement 


11 


provided temperature control of the sealed fluid. A transparent housing 
permitted visual observations on the rate of gas flow out of the pressurized 
cavity; and, in studies on grooved rotors, a stroboscope permitted obser- 
vations of the gas ingestion process in the rotating grooves. 

Figure 4 is a diagrammatic sketch of the experimental apparatus used 
for evaluating viscoseals in sodium . The temperature was maintained by 
a furnace (resistance heaters) surrounding the viscoseal housing. The 
temperature was monitored by thermocouples attached to the seal housing 
outer diameter and by one thermocouple submerged in the sodium in the 
pressurized cavity. 

The enclosure surrounding the viscoseal assembly is pressurized 
with argon (after prior evacuation) to slightly above ambient to insure 
exclusion of air. Prior to the introduction of liquid sodium, the input 
shaft is set at some fixed speed. Sodium is introduced by pressurizing 
the reservoir and venting the leak detector tank (see sodium supply system 
in fig. 5). The sodium first fills the pressurized cavity (viscoseal high- 
pressure end) and then is allowed to fill the leak detector tank to a pre- 
determined level. The leak detector tank float displacement is monitored 
by recording the differential transformer output that is produced by the 
transformer core attached to the top of the float. Thermocouples placed 
at different levels in the tank provide calibration during the filling process. 

By holding the reservoir temperature at 220° F, a low degree of oxide 
solubility was maintained. The 20-micron filter (stainless steel) inside the 
tank allowed 220° F filtration. A second filter external to the reservoir 
(5 microns) provides additional filtering at approximately 300° F. Sodium 
is partially removed from the system by pressurizing the leak detector and 
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forcing the sodium back into the reservoir. By alternately filling the 
system and then partially returning the sodium to the reservoir, the so- 
dium can be recirculated and refiltered. This recirculation provides a 
hot flush or cleaning action on the viscoseal assembly before a run. (All 
transfer lines and the test section are held at 300° to 500° F during this 
recirculation process. ) For all sodium studies the leak detector tank was 
held at 300° F, and any seal leakage resulted in a drop of liquid level and 
float position. Gas ingestion into the pressurized cavity was indicated by 
a rise in sodium level in the leak detector. 

Sealing capacity, when sealing sodium, was obtained by increasing 
cavity pressure until the leak detector indicated leakage. Line thermo- 
couple readings were also found to be sensitive leak indicators. Visco- 
seal geometries employed in this evaluation are shown in table H. The 
geometry near the optimum (derived in ref. 19) was taken as the basis 
for comparison, and is listed in the first row of table n as geometry 1. 
Secondary grooves were added to the lands of this basic geometry to form 
the second geometry evaluated. For the third geometry, a ratio of groove 
to land width of 1. 90 was selected; and for the fourth geometry, secondary 
grooves were added to geometry 3. Geometry 5 is the internally grooved 
housing having the same groove dimensions as geometry 1. The rotor 
outer diameters were 1. 993 inches, and the housing bores were 2. 001 
inches. 

EXPERIMENTAL RESULTS AND DISCUSSION 
Pressure Patterns 

Figure 6 shows an experimentally determined pressure pattern in a 
viscoseal composed of an internal helically grooved housing and a smooth 
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rotor arrangement. In a plane orthogonal to the rotor centerline, the 
pressure increases across the groove (from point A to B) in the direction 
of rotor rotation and decreases across the land (from point B to C). The 
increase and decrease in pressure repeat for each groove-land pair and 
produce a saw-tooth pressure profile in the orthogonal plane. The pres- 
sure increases from the low-pressure end (plane 3) to the high-pressure 
end (plane 1), and the pressure pattern has a helical twist corresponding 
to the helical grooves. (For a rotating grooved shaft this pressure pattern 
would be rotating. ) 

The pressure gradients in the planes orthogonal to the rotor center- 
line are essentially linear for axial positions not influenced by the end 
effects. Figure 7 shows typical experimental results. These pressure 
profiles are the same as those for the orthogonal planes shown in figure 6 
except that a single groove -land pair is unwrapped to form a plane figure. 
(Neglecting curvature is not a significant error since the ratio of rotor 
radius to groove depth is 65 to 1. ) 

In these orthogonal planes, the ratio of pressure to cavity pressure 
P/P c increases linearly over the groove and then decreases linearly over 
the land. This linear relation was found to hold to the maximum test speed 
of 6000 rpm, which produced a modified Reynolds number (Re*) of 0. 25. 

End Effect 

When the helical grooves connect directly with the pressurized 
cavity, the pressure developed along the land leading edge decays sharply 
near the high-pressure end of the seal. Typical results are shown in fig- 
ure 8 for 1000 and 6000 rpm. The pressure pattern at the end was essen- 
tially the same for the speed range investigated (1000 to 6000 rpm) except 
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that the results for 5000 rpm (not shown) and 6000 rpm start to show a 
slight shortening of the end-effect length. This shortening is probably due 
to convective inertia effects. In all cases, the end-effect axial length was 
approximately equal to one-half the axial groove width. This suggests that 
end-effect length could be taken to one-half groove width without significant 
error. The pressure decay within the groove is due to the equalization of 
the land leading- and trailing-edge pressure, since a common pressure 
must exist in the cavity. 

Gas Ingestion and Sealing Capacity Observations 
When Sealing Water 

Grooved housing . - The series of photographs in figure 9 shows the 
gas ingestion process when sealing water with a grooved housing and 
smooth outer diameter rotor. At a modified Reynolds number (Re*) of 
1. 75 (2000 rpm), no gas ingestion is evident, the water film is clear, the 
interface is stable, and the scavenging length is nonwetted. 

At Re* = 2. 58 (3000 rpm) the film had large gas pockets, gas was 
ingested as evidenced by air bubbles passing out of the transparent vent 
line, and the sealing capacity (parameter, XG) was lower than that at 
Re* = 1. 75 (2000 rpm). Various degrees of gas ingestion and film rupture 
are evident for Re* = 3. 50 (4000 rpm) to 6. 10 (7000 rpm). At Re* = 7.00 
(8000 rpm), no detectable gas ingestion rate was evident, and gas ingestion 
was not detectable from Re* = 7. 00 (8000 rpm) to the highest speed evalu- 
ated, which gave Re* = 10. 84 (12 000 rpm). At Re* = 7. 00 and above, 
the seal wetted length had two distinct regions (fig. 10(f)): a gas-liquid 
region, and a region of homogeneous liquid at the high-pressure end. 
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Between Re* = 3. 5 and 6. 10, it was observed that the bubble size 
progressively decreased when moving from the low-pressure end to the 
high-pressure end. It is probably that with longer seal lengths the corre- 
sponding higher pressures would reduce the bubble size sufficiently to al- 
low escape over the lands, and thus gas ingestion would be eliminated. It 
was observed that concentricity of rotor and housing markedly affect gas 
ingestion occurrence. 

Grooved rotor . - In a series of visual observations using a grooved 
rotor, the gas ingestion started at Re* = 1.9 and increased in rate with 
increasing Reynolds number (maximum investigated, Re* =10.84). Below 
Re* =1.9 the liquid-to-gas interface was stable and no gas ingestion was 
observed. 

Gas ingestion mechanism . - These observed gas ingestion character- 
istics of the grooved rotor and the grooved housing are illustrated in fig- 
ure 10. The probable mechanism producing the difference between the gas 
ingestion process of the grooved rotor as compared with that of the grooved 
housing is illustrated in figure 11. For the grooved housing, when the cen- 
trifuge action on the liquid becomes great enough, the gas bubbles are dis- 
placed to the rotor surface and pass over the lands. Therefore, the bub- 
bles are not pumped to the high-pressure end. For the grooved rotor, the 
centrifugal action on the liquid forces the gas bubbles into the grooves; 
therefore, the bubbles are pumped to the high-pressure end of the seal. 

When a viscoseal is ingesting gas, the axial pressure gradient is non- 
linear and increases as the liquid-to-gas volume ratio increases in the di- 
rection toward the high-pressure end. For no gas ingestion, the axial 
pressure gradients are linear as shown in figure 8, which contains data 
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on sealing oil. As shown in figure 12, gas ingestion results in nonlinear 
pressure gradients (increasing slope with increasing liquid-gas volume), 
which are significantly different from the gradients shown in figure 8. 

Leak rate. - In both viscoseal combinations (grooved housing and 
grooved rotor) water was sealed, with no detectable leakage. It was ob- 
served that, although gas ingestion reduced the sealing capacity AG, ef- 
fective sealing of the water was maintained. 

Results of Sodium Sealing Studies 

Grooved housing . - In sealing sodium, the grooved housing did not 
ingest gas (fig. 13) for concentric alinement (0. 005 in. F. I. R. ) over the 
range investigated (Re* = 8. 66 to 50. 5) and for pressures ranging from 
2 to 100 psig. The sensitivity of the leak detector to measure gas inges- 
tion rate was estimated to be 0. 50 cubic centimeter per hour. It was ob- 
served that nonconcentricity leads to gas ingestion, but this situation was 
not investigated. 

Grooved rotor . - The grooved rotor started to ingest gas at Re* = 10, 
and the rate increased with increasing Reynolds number, as shown in fig- 
ure 13. This ingestion characteristic is similar to that observed when 
sealing water with the grooved rotor. 

Gas ingestion into closed cavity . - In sealing sodium, gas ingestion 
was readily observed by monitoring the liquid level in the leak detector. 

As shown in figure 14, a straight line on the oscillograph paper indicated 
no gas ingestion, no leakage, and thermal equilibrium. When gas inges- 
tion occurred, the liquid level line showed an upward trend, the slope of 
which is the ingestion rate. As gas ingestion proceeds, the gas bubble in 
the pressurized cavity increases in size (fig. 14(a)), and when the size 
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increases to the rotor outer diameter, the gas flows back out the annulus 
formed by the rotor and housing. This causes a decrease in bubble size 
and a sharp drop in the leak detector readout trace. If the interface is far 
enough away from the seal low-pressure end, the liquid will be scavenged 
back within the nonwetted length, and no loss of liquid will occur. The in- 
gestion process will begin again and the net result is a saw-toothed trace, 
each sharp drop indicating gas blowback. If the interface is sufficiently 
close to the seal end, liquid loss will occur, and the leak detector readout 
will show a saw-toothed trace with a downward trend, which represents 
a liquid loss at each blowout (fig. 14(c)). This latter action may be the 
’’seal breakdown” phenomenon reported in reference 21 and the "secondary 
leakage” phenomenon reported in reference 4. (In sealing water, the re- 
peating bubble growth and subsequent seal blowout was visually observed 
by means of a transparent housing. ) 

Sealing capacity . - The grooved housing had a higher sealing capacity 
(higher sealing parameter, XG) than the grooved rotor. The comparison 
is given in figure 15, which gives the sealing parameter XG as a function 
of Reynolds number based on c (Re ) . The difference between the value of 
XG for the grooved housing and the value for the grooved rotor was at- 
tributed to the difference in gas ingestion characteristics. The grooved 
rotor operates with a gas-liquid mixture along the full wetted length simi- 
lar to that shown in figures 9 (b) to (e); the grooved housing operates with 
a homogeneous liquid at the high-pressure end of the seal and a gas-liquid 
mixture at the low-pressure end similar to that shown in figure 9(f). 

Figures 16(a) and (b) show the sealing parameter XG as a function of 
Reynolds number (Re ) for groove configurations of five- and ten-helix 
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starts (see table II) with and without secondary grooves on the lands. In 
both cases, a higher sealing coefficient was obtained through the use of 
secondary grooves. Secondary grooves also operate at lower power ab- 
sorption levels for equal sealing capacity, as was determined by the 
amount of heating power required to maintain operating temperature. 

This lower power absorption is due to less land area as compared with 
the land area of the seals without secondary grooves. 

Sodium was sealed for pressures ranging from 2 to 100 psig with 
negligible leakage loss. The operating range was within Reynolds numbers 
(Re*) of 8. 66 to 50. 5, and the fluid temperature was in the range 300° 
to 625° F. The maximum length of any one run was 8 hours and accumu- 
lated testing time amounted to 50 hours for 10 runs. The oxygen content 
of the sodium for the 10 runs varied between 50 to 95 parts per million, 
as determined by the analytical method described in reference 31. 

THEORETICAL ANALYSIS OF THE VISCOSEAL 
Basic Model 

The rectilinear Cartesian coordinate system was selected for mathe- 
matical tractability. A parallel groove-ridge plate geometry with a 
smooth flat plate moving relative to the parallel groove -ridge plate, both 
of infinite horizontal extent, was chosen as the basic model (see fig. 17). 

Conceptually this model can be thought of as an infinitely long smooth 
rotor with a very large diameter and a concentric grooved housing. Actu- 
ally this infinite grooved plate model is the limiting case of a helical 
grooved cylinder where R — 00 or c — 0. If c/R<0.01, this model 
should give a good qualitative picture of the flow field and pressure gener- 
ation and significant trends for optimum geometry. (See argument in 
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appendix A. ) For example, a practical geometry is a 2-inch-diameter 
rotor with a 0. 005 -inch concentric clearance, c/R = 0. 005. 

The clearance between the infinite plates is filled with any incompres- 
sible fluid, and it should be noted that for this infinite horizontal extent 
case there are no end effects (in practical engineering design, of course, 
end effects must be accounted for). Since experiments at Lewis show that 
the end effect zone decreases with increasing rotor speed and also that it 
is less than a groove width, this analysis should be valid in the bulk of the 
homogeneous fluid (see figs. 18 and 19). 

In the flat -plate case, the boundary conditions on a groove -ridge plate 
in motion are constant; however, it will be advantageous to have the 
smooth flat plate move relative to the grooved flat plate. This will elimi- 
nate the temporal inertia terms and reduce the subsequent independent 
variables by one. 

A word of caution for extension of the results to a case where the 
groove -ridge pairs on a finite shaft are rotating relative to a fixed smooth 
housing. The boundary conditions are different when the smooth housing 
is stationary and the grooved rotor is moving, and when the housing is 
grooved and the rotating shaft is smooth. When the rotor is grooved, the 
surface speed is different at each point along the step or edge of the moving 
boundary. Thus this analysis should be more valid for a smooth rotor and 
grooved fixed housing (from ref. 5 the grooved housing appears to be the 
desired mode of operation) with a large radius of curvature and/ or small 
clearance. 

When referring to figure 17, it is noted that by the choice of coordinate 
system the groove-ridge pairs appear only in the y, z-plane. Later it will 
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be shown that this choice of coordinate system will be necessary to reduce 
the analysis to a two-dimensional problem for the quasi two-dimensional 
flow field. It should also be noted that the drag force has been resolved 
into its components in this coordinate system and there is no equilibrium 
of forces in the z-direction with inertia (nonlinear) effects neglected. 

BASIC EQUATIONS 


The Navier-Stokes equations for a homogeneous, incompressible, 
laminar, Newtonian fluid are the following (see ref. 32): 
x-direction: 
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+ J*.* + A + 

& Vax 2 ay 2 0z 2 . 


(i) 


y-direction: 




/_2 _2 -.2 ' 
/a_v + ^v + a_v 

Vax 2 3y 2 0z 2 


( 2 ) 


z-direction: 


P 


Dw 

Dt 




The Eulerian or spatial derivative is 
local acceleration (temporal velocity 
eration (spatial velocity change): 


a^ + a^ + AA 

.ax 2 ay 2 az 2 / 

D/Dt, which is composed of the 
change) and the convective accel 
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Dt 




d_ 

at 


a a a 

+ U + V + W 


■J V. 


dx 


dy 


dz 


J 


Local, unsteady, or 
nonstationary accel- 
eration 


Convective 

acceleration 


me incompressible continuity equation is 

ii! + iz + ^ = o (4) 

3x 3y dz 

The incompressible energy equation with constants C and k and neglect- 
ing compression work DP/Dt (ref. 32) 

pgc5I = k^ + A + 3V| +jw + Q . 


Dt 


\3x^ dy 
where <& is the dissipation function 


(5) 


az 4 


= 2 ffinf + (M 2 + + (?i + isf + 

l\ax/ Vay/ Vaz/ J Vax ay/ V<ty az/ 


/ 3u + 3w\^ _ 2 /3u + 0v + 3w V 
Vaz 3x / 3 \3x dy dz / 


The assumptions for this analysis are as follows: 

(1) For steady flow, 


3u _ 3v _ aw _ q 

at at at 


(2) For no body forces, 


F = F = F =0 
x y z 


This means there are no electromagnetic fields present and gravitational 
effects are negligible. 

(3) Isothermal and equal temperature plate surface conditions are 


assumed. 
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(4) Viscous (frictional) heating is negligible. Therefore, the thermo- 
physical properties are constant. This appears to be a good assumption 
for order 10 (water) and smaller (liquid metals) Prandtl number Pr fluids: 



k 


hi experimental setups there is an inherent unsteadiness due to vibra- 
tions, shaft runout, and eccentricity of rotor with respect to housing and 
misalinement. These factors will influence the degree of correlation be- 
tween the analytical and experimental results. 

The conservation of momentum and mass equations were nondimen- 
sionalized and a formal ordering procedure was used to determine the 
relative magnitude of the terms. Consequently, in a formal way the neg- 
ligible terms were found (see appendix A). 

An important parameter, the modified or reduced Reynolds number, 
was found: 

2 

Re * =Re L(f;) =Ee c(i;) * C/L<1 
This modified Reynolds number gives the relative magnitude of the inertia 
forces to the viscous forces; that is, 

Re* = forces 

Viscous forces 

Thus, inertia forces can be neglected only if Re* < 0. 01. 

After the formal ordering procedure (see appendix A), the three- 
dimensional flow field equations are as follows: 
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Conservation of momentum: 


0U 0U 0U 1 0P 0^U 

U + V — + W — = + V 

0x 0y 0z p 0x gy2 

(x-direction) 

(6) 

— = 0 (y -direction) 

3y 


(7) 

u5w + v SE + w iE = _ l 8P | | aV 

ax By 3z p 3z \g y 2 az 2 ; 

j (z -direction) 

(8) 

Conservation of mass: 



0u + 0v + 0w _ o 
0x 0y 0z 


(4) 


Fully developed flow is assumed in the x-direction. This means that 
V and 8P/9x are not functions of x and that the pressure gradient in the 
x-direction (parallel to the groove) is a constant (verified experimentally 
in ref. 5). In figure 18, it is observed that this assumption is valid only 
where the average pressure profile at the ridge leading or trailing edges 
in the axial direction is linear. Gas ingestion results in a nonlinear pro- 
file (see figs. 18 and 19). 

With the assumption of fully developed flow along the groove-ridge 
direction (x-direction), the flow field equations become two-dimensional: 

v — + w — = -- — + v iLH where — = constant (9) 

dy 9z p 0x 8 y2 0x 
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Se! 0 

3y 

(10) 


= -i *Z + v(&L + i2\ 

(ID 

3y 9z 

p 0Z \0y 2 0z 2 / 

Conservation of mass: 


^ + ^ = 0 
3y 3z 

(12) 


The boundary conditions for an arbitrary shaped groove -ridge geome- 
try (see fig. 20) are: 

( 1 ) 

v = 0 ''I 

u = U cos a > at y = 0 
w = -U sin a J 

(2) 

u = v = w = 0 at y = h(z) 

For example, when h(z) is a rectangular groove-ridge pair the boundary 
condition is: 

u = v = w= 0 at c < y < c + h Q (on side of groove) 
u = v = w = 0 at y = c (on ridge) 
u = v = w = 0 at y = c + h Q (on groove root) 

(3) Periodicity of groove -ridge pairs at the clearance interface. 

For 0 < y < h(z) 

v(0, y) = v(L sin a, y) 
w(0, y) = w(L sin a, y) 

P(0, y) = P(L sin a, y) - (AP r+g ) z 
where (AP r+g ) z = const 
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All velocity derivatives must also be the same to satisfy the periodicity 
boundary conditions 

a n+m v(0,y) _ a n+m v(L sin a, y) 

0y n 3z m ay 11 3z m 

3 n+m w(0, y) _ 8 Ili ~ lu w(L sin a, y) 
ay 11 az m ay 11 az m 

for 0 < n, m < » 

(4) Boundary condition for aP/3x for a cylindrical geometry (visco- 
seal). 

In appendix B it is shown that for a cylindrical geometry 3P/3x de- 
pends on the (AP r+ g) obtained in the z -direction 

3P _ N(AP r tg >z 
3x 27 tR cos a 

These equations (9), (11), and (12) can be solved for the velocity com- 
ponents. The velocity components can be obtained for arbitrary (AP ) . 

I+g It 

Then the axial pressure gradient 3P/ 3L’ can be found from 

3P _ AP axial _ ( N )( AP r+g ) z CQt a 
3L’ L’ L 

where L' is the axial length (see fig. 21). Qualitative pressure distribu- 
tion from experimental data in the x- and z -directions is shown in figures 
22 and 23. 

By examining the analysis, boundary conditions, and the derived equa- 
tions, the axial pressure gradient aP/3L’ is a function of (a, U, p, v), 
where cr is a shape parameter that is a function of (b/h Q , c/h Q , a/b, a). 
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Since convective inertia effects are important, new geometric parame- 
ters appear which are not implied in creeping flow theory. Most notable of 
these is the aspect ratio b/h Q , which implies the degree of convective ef- 
fect. When the b/h Q and a/b ratios are known, N, the number of helix 
or groove starts, can be found. (Direct effect of N is seen from boundary 
condition (4). ) Thus the number of groove or helix starts also is an impor- 
tant parameter. Since inertia is important, a high density is desired for 
pressure generation. 

The equations (9) to (12) and boundary conditions are the minimum 
number of terms that are necessary to mathematically describe the physics 
of the internal flow. This is the simplified form of the Navier Stokes equa- 
tions which should be used for an analytical solution; however, a mathe- 
matical function to describe the groove -ridge shape may be too complex 
for a closed-form analytical solution. Thus the equations are solvable by 
using high speed digital computer numerical methods. Since numerical 
schemes can accommodate all two independent variable terms with equal 
ease, it is suggested that the complete set of two-dimensional equations be 
solved numerically, i. e. , 




By 


3w 

v — + w 


3z 

0W 


3u 3u 

v — + w — 

3y 3z 


3y 3z 

1 3P 


p By 

1 3P 

p 3z 


\.3y 3z 
+ + s2w ' 


\dy‘ 


dz‘ 


p 3x 


i 2 S 
Z 3z 


3P 


where — = constant 
3x 


(13) 


(ID 


(14) 


^ + ^=0 
3y 3z 


( 12 ) 
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Of course, the numerical solution will determine under what condi- 
tions the previously derived simplified form of the two dimensional Navier 
Stokes equations are good approximations. 

By examining the set of equations (11) to (14) it is seen that v and w 
can be calculated from equations (11), (12), and (13). The solution is 
identical to a two-dimensional flow case. Then the u flow field can be 
found from equation (14) which is linear in u. This is an example of the 
Independence Principle. 

From appendix A the energy equation with no temperature variation 
along the groove-ridge direction (3T/3x = 0) is 



d8 _ v j 

f a 2 e + 

aV\ 

| + JL 


(15) 

Sy 

dz Pr ' 

w 2 

8z 2 ^ 

1 PgC 

pg c 


Comparing the above equation with the x-direction momentum equation 
it is seen that the mathematical form of both equations is identical. 
Substituting 6 =^u 

and 

— — $ or — — = a constant =2> - - — 

PgC pgC p 3x 

Thus with Pr = 1 and the above substitutions, the equations (14) 
and (15) are identical. The thermal boundary conditions are 
At y = 0 

8 = constant =^>w --U sin a = constant 

At y = h(z) 

6 = 0=>w = 0 

and the boundary conditions are similar in form. 

The above statements mean that once the u velocity field is found 
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the solution will also give the temperature distribution for the stated re- 
strictions of constant frictional heating or a heat source in the fluid with 
the appropriate physical parameters equivalent to the constant pressure 
gradient along the groove-ridge pair direction and Pr = 1. Note that this 
temperature distribution solution violates the initial assumption of con- 
stant thermophysical properties; however, for small temperature differ- 
ences the above solution will give a good qualitative picture of thermal 
effects. 

DISCUSSION 
Discussion of Equations 

This analysis is valid only for a homogeneous fluid and figures 18 
and 19 clarify the region of validity. The analysis is not valid for a 
scavenging area and a region consisting of a mixture of gas and liquid 
(see fig. 18) or under conditions of gas ingestion. 

It should also be pointed out that the analysis is not valid in the end 
effect region at the high pressure end of the seal (see fig. 18) where the 
groove connects with the cavity. Since the appropriate equations are non- 
linear, the pressure is not a harmonic function; thus, the usual analog 
methods of finding end effect are not available as found in creeping flow 
regions (see appendix C). 

The equations are unique in many ways. Mathematically, the problem 
has been reduced to two dimensions (y- and z -directions) but all three ve- 
locity components are present. Thus the flow field can be described as 
quasi-two-dimensional. Also, there is a pressure gradient in the third 
dimension (x -direction), but fortunately this gradient is a constant and 
dependent on the y- and z -direction gradients. 
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The velocity component v is no longer negligible as it was in the 
creeping flow solution. Physically, v cannot be negligible because of the 
presence of the edge effect which will result in a component of velocity in 
the y-direction. 

Since there is no restriction on net volume flow, the pressure gen- 
eration prediction equations should equally apply to hydrodynamic seals, 
bearings, and pumps. 

When further examining the physics of the problem, it is seen that 
across the groove-ridge pairs the unbalance of the convective accelera- 
tion force with the drag force results in a pressure gradient in the 
z-direction. Since the step or edge causes the convective acceleration, 
its effects are not negligible unless there is a very large aspect ratio 
and the mathematical model accounts for the edge effect. Along the 
groove and ridge (x -direction) the drag force is in equilibrium with the 
induced pressure and convective inertia force. The drag force not only 
causes fluid motion but also aids in maintaining equilibrium with the 
pressure along the groove or ridge. 

Both equations of motion illustrate the coupling of u, v, and w in 
the convective inertia terms. This means that physically a condition such 
as no end leakage is possible. In creeping flow there is no coupling be- 
tween the velocity components, which illustrates the shortcoming of that 
analysis (i. e. , end leakage is not zero). 

The assumption that the inertia terms are of the same order as the 
viscous terms results in equations of the boundary layer type; however, 
the boundary conditions are different. 
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The equations are of a form that might tempt one to reduce the num- 
ber of independent variables by one by trying to find the admissable flow 
field which would result in a similarity transform solution. This tech- 
nique fails because of the characteristic lengths that are present in this 
type of confined flow. 

The equations are solvable by using high-speed digital computer nu- 
merical methods. A mathematical function to describe the groove-ridge 
shape may be too complex for a closed-form analytical solution. It 
should be interesting to find out if the step is still the optimum shape 
with convective inertia effects included. For creeping flow, Lord Rayleigh 
found in 1912 that the step was the optimum shape by using calculus of vari- 
ations. 

The equations and ordering procedure should be reexamined carefully 
before extending these results beyond the assumptions used in this analysis. 

Discussion of Modified Reynolds Number 

A check was made on the significance of the modified or reduced 
Reynolds number by calculating Re* from the available published data at 
the claimed onset of turbulence (breakpoint). In table El it can be noted 
that the departure in all cases from a constant sealing coefficient occurred 
when Re* was in the range of 0. 26 to 2. 2, which suggests that the so- 
called beginning of turbulent operation is really the beginning of significant 
convective inertia effects. This would be in agreement with Kettleborough 
(ref. 29), who found qualitative agreement between his inertia solution and 
experimental turbulent results. 

For seals 5, 6, and 7 (see ref. 27 and table El), only the groove to 
ridge ratio has changed. As the ratio increases, the breakpoint decreases. 
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This is reasonable since a larger ridge area results in a less influential 
convective inertia effect and delays the beginning of pronounced convective 
inertia effects. It can also be seen in table III that the calculated Re* are 
within an order of magnitude of one another for all of the experimenters. 

It is difficult, nevertheless, to select the breakpoint (see fig. 24). Even in 
a range of possible breakpoint numbers, however, the order of magnitude 
will remain the same. Further experimental results will be necessary, but 
in general, it can be concluded that if Re* > 1, a constant sealing coeffi- 
cient or bearing number will not exist. 

In figure 25 it is seen that classically the empirical sealing parameter 
XG was thought to be a constant until some critical Re and then to vary 
exponentially in the turbulent regime. The creeping flow solution was 
thought to be valid until the critical Re and then a semiempirical turbulent 

V 

theory had to be used. Actually, the onset of significant convective inertia 
may be gradual, as shown in figure 25. Data (from ref. 27) for seal 7 in 
figure 24 tends to show this. The derived equations herein pertain to all 
cases under assumptions stated and will break down only under strong turbu- 
lence or an unknown instability phenomena. However, even in the turbulent 
flow model the convective inertia terms must be considered in addition to 
the Reynolds' stress terms. This is specifically shown by Kulinski and 
Ostrach (ref. 33) in their critical evaluation of high speed fluid film lubrica- 
tion theory. 

Principle of Viscoseal Operation 

The classical explanation of viscoseal operation is now described. 

Since the seal is "pumping" the sealed fluid as fast it is leaking out, there 
is no net end flow (zero leakage). From the basic laws of fluid mechanics, 
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however, it is seen that the viscoseal axial pressure gradient results from 
the unbalance of the convective inertia force across the groove -ridge pairs 
and the drag force. The convective inertia force is much greater than the 
opposing drag force, and this results in a net pressure force. Along the 
groove -ridge direction, the drag force maintains equilibrium with the pres- 
sure force and convective inertia force. All three velocity components are 
"coupled, " and thus zero net leakage is achieved. This is an excellent ex- 
ample of the physical importance of a nonlinear effect. The inadequacy of 
the linear theory (creeping flow) has been previously discussed. 

SUMMARY OF RESULTS 
Pressure Patterns 

Pressure patterns were studied for a viscoseal in the laminar-flow 
regime with a mineral oil and the following results were obtained: 

1. For the condition of helical grooves in direct connection with the 
pressurized cavity, the end effect is evidenced by a sharp decay of land 
leading-edge pressure. 

2. No significant change in end-effect length was noted over the speed 
range investigated (modified Reynolds numbers (Re*) of 0.04 to 0. 25). 

3. Circumferential pressure gradients in the plane of rotation (plane 
orthogonal to 77 -axis) increase linearly across the groove and decrease 
linearly across the land. This pressure pattern repeats for each groove- 
land pair and provides a saw-tooth pressure profile around the circumfer- 
ence. The pressure gradients along the groove and land were linear. 

Gas Ingestion and Sealing Capacity 

Visual observations and experimental evaluation of viscoseal assem- 
blies with water (modified Reynolds numbers of 1. 75 to 10. 84) and sodium 
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(modified Reynolds numbers of 8. 66 to 50. 50) as sealed fluids disclosed 
the following: 

1. Gas ingestion characteristics were fundamentally different for the 
helically grooved rotor-smooth bore housing combination. In the grooved 
housing- smooth rotor combination, centrifugal action tended to hold the 
gas bubbles out of the groove; thus the bubbles tended to escape over the 
lands rather than be pumped into the seal high-pressure end and a condi- 
tion of no gas ingestion could be achieved. In the grooved rotor -smooth 
bore housing combination, centrifugal action forced the gas bubbles into 
the rotor grooves; thus the gas bubbles were pumped into the seal high- 
pressure end and increasing speeds resulted in increasing gas ingestion. 

2. In sealing sodium, the grooved housing- smooth rotor combination 
had no detectable gas ingestion rate, but the grooved rotor-smooth bore 
housing combination showed increasing gas ingestion with increasing 
Reynolds number. 

3. No gas ingestion occurred when sealing water with either the 
grooved rotor-smoothbore housing or the grooved housing- smooth rotor 
viscoseal when the liquid interface was stable (below Reynolds number (Re*) 
of 1. 9). 

4. In sealing sodium, the grooved housing- smooth rotor combination 
produced a higher sealing capacity (higher sealing parameter) than the 
grooved rotor-smooth bore housing combination. This was attributed to 
the presence of gas ingestion in the case of the grooved rotor. 

5. Sealing, with negligible sodium leakage rate, was obtained for 
short operational periods (maximum length of continuous evaluation was 
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8 hr) over a pressure range of 20 to 100 lb force/ sq in. gage and at sodium 
temperatures of 300° to 625° F. 

Theoretical Analysis of the Viscoseal 
A review of the literature on the theoretical analysis of the viscoseal 
revealed that either the convective inertia forces or the viscous forces 
were neglected (inviscid fluid assumption) . This analysis included both 
convective inertia and viscous forces. The following two-dimensional set 
of equations for a quasi-two-dimensional flow field were formulated and 
are solvable for the velocity and pressure distribution and optimum geom- 
etry utilizing numerical methods on a high-speed digital computer: 
Conservation of momentum: 

v — + w — = - i — + v where — = constant (9) 

3y Bz p dz gy 2 9x 

-s0 (10) 

3y 


aw 


+ w 


aw 


ay az 


1 3P 

p 3z 
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Conservation of mass: 


+ = 0 
By Bz 

With the following boundary conditions 

(1) At the moving smooth plate surface, 

u = U cos a v = 0 w = -v sin a 

(2) At the fixed parallel groove geometry surface, 


(12) 


u = v = w = 0 
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(3) Periodicity of groove-ridge pairs at the clearance interface. 
For 0 < y < h(z) 


v(0, y) = v(L sin a, y) and 


a n+m v ( 0 , y ) _ a^'nytL sin a,y) 
ay” 9z m 9y“ 0z m 


nXTTI v ^n + ITl • \ 

w(0, y) = w(L sin a, y) and - 1 w<,L sin a, y; 

ay 11 0z m ay 11 az m 

P(0, y) = P(L sin a, y) - (AP r+g ) z for 0 < n, m < «> 

( 4 ) 

ap . N < AP r + g>z 

3x 27rR COS Of 

From the development and analysis of the equations, the following 
results and conclusions can be made: 

1. A modified or reduced Reynolds number evolved which gives the 
relative magnitude of the convective inertia forces to the viscous forces: 


Re* = Re^ 




if 


* <1 

L 


The characteristic length was selected to be the length of the groove - 
ridge pair in the direction of the smooth plate velocity vector. This is 
the principal length over which the driving force acts. The flow behavior 
over each groove -ridge pair is the same. Convective inertia effects can 
only be neglected if Re* « 1, which will result in true creeping flow. 

2. A modified Reynolds number equal to 1 will predict a significant 
change in sealing coefficient or bearing number which heretofore was at- 
tributed to the onset of turbulence. This means that the convective inertia 
effects are the foremost contributing factor to the improvement in sealing 


i 
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coefficient or bearing number, not the turbulence effects. A calculation 
of the modified Reynolds number at the onset of turbulence of four different 
experimenters disclosed that the modified Reynolds number was near 1 
(Re* = 0. 25 to 2. 2). Thus a turbulent flow analysis must also include the 
convective inertia terms in addition to the appropriate Reynolds stress 
terms. 

3. Several additional parameters of pressure generation dependence 
were found that are not implied in the creeping flow analysis. These addi- 
tional parameters, which arise from convective inertia effects, are the 
fluid density p, the groove width to depth or aspect ratio b/h Q , and the 
number of groove-ridge pairs N. 

4. Solving equations (11) to (14) by numerical analysis will permit the 
investigation of arbitrary groove-ridge shape, for example, sinusoidal, 
rectangular, triangular, etc. Now a check can be made to see if a step is 
the optimum shape with convective inertia effects as it is for creeping flow. 

5. The analysis provides a unified approach to parallel groove opera- 
tion previously covered by creeping flow (sometimes called laminar flow) 
analysis and a semiempirical turbulent flow analysis. 

6. No restriction was placed on the net flow situation between the 
plates. Therefore, the analysis can be used to analyze parallel groove 
geometries in hydrodynamic seals, bearings, and pumps. 

7. The viscoseal axial pressure gradient results from the unbalance 
of the convective inertia force across the groove-ridge pairs and the drag 
force. Along the grooves or ridges, the drag force maintains equilibrium 
with the pressure force and the convective inertia force. 
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8. Zero leakage in viscoseal analysis can be attained through the con- 
sideration of convective acceleration terms (nonlinear effect). It is argued 
that the creeping flow analysis, which is a linear theory, cannot prove zero 
leakage. 

9. The format of this analysis can be used to evaluate qualitative con- 
vective inertia effects in any step or converging wall hydrodynamic seal, 
pump, or bearing (Rayleigh step, journal bearing, etc.). 

10. Solution of the u velocity (along the groove -ridge direction) will 
give an insight to the heat transfer in the viscoseal for the case when 
Pr = 1. 

APPENDIX A 

FORMAL ORDERING PROCEDURE FOR SIMPLIFYING 

BASIC EQUATIONS 

The Navier -Stokes equations for a homogeneous incompressible 
Newtonian fluid, neglecting body forces and for steady flow, are 



The incompressible continuity equation is 

^ + l}! + ^ = 0 (A4) 

8x 9y 9z 

To find the important terms in the previous flow field equations, a 
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formal ordering procedure will be used to determine the relative magni- 
tudes of the terms. 

The terms are normalized by scaling them to their appropriate maxi- 
mum characteristic value as follows: let 

x = — x - y = Z z = 

L cos a c L sin a 


U cos a V U sin a 

The characteristic dimensions c and L can be seen in figure 17. 

It is noted that a trigonometric relation exists between x and z charac- 
teristic lengths and that c/L < 1. The characteristic length was selected 
to be the length of the groove-ridge pair in the direction of the smooth 
plate velocity vector. This is the principal length over which the driving 
force acts. The flow behavior over each groove -ridge pair is the same. 

Let p* and v* be the fluid density and viscosity corresponding to a 
reference state (P*, T*). Therefore, 

V = — = 1 

l/* 


and 

p* 

for incompressible fluids with constant thermophysical properties. Now 
to find the relative magnitude of V so that all terms of the incompressible 
continuity equation will be of the same order, the incompressible continuity 
equation is nondimensionalized: 


3u + V L 3v + 9w 
3x U c 9y Bz 


(A5) 
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To have all terms of the equation the same order of magnitude, ^ must 

be of the order of 1. Therefore, V = U — , and thus "v = — — . 

L U c v 

Next the equations of motion are nondimensionalized. The pressure 
term must first be nondimensionalized; even though the pressure is not 
the driving force, it is the function of primary importance. Thus the char- 
acteristic pressure is not extremely important. Since the pressure gradient 

is related to the rotor speed, it appears to be logical to nondimensionalize 

o 

the pressure with respect to the dynamic pressure pU (also called the 
dynamic head): 



2 

It can be seen that the characteristic dynamic pressure pU has re- 
duced the number of nondimensional groups by one (i. e. , if P* (a refer- 
ence pressure) had been used, another nondimensional group would be 
present in the equations). Also, since the incompressible case is being 
analyzed, the pressure appears in the dynamic role only, not in the thermo- 
dynamic role. 

Nondimensionalization of x-Direction Momentum Equation 


When the x -direction momentum equation is nondimensionalized, 
equation (1) results in 
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Next a Reynolds number is defined based on the characteristic length L: 

UL 


Re L " 


v* 


where 


and thus 


— 3u 9u — 9u 
u — + v — + w — 


v = v*v = v* 


0P 


0X 


3y 
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0x 

cos a 
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,2 ^ 2 - 


cos a 0x 


u( ^yo + _i 

2 Vc J a _2 


0 2 u 


2 —2 
sin a 0z 


(A6) 


The convective acceleration and pressure terms are of unit order; thus, 

the viscous term must be of unit order also. The largest term of the 

2 2 - 2 
viscous forces is (— \ since f— ^ < 1; therefore, the — i ^ 

\c ) 0 -2 \lJ Re L \c ) 


should be of unit order. The modified Reynolds number is defined by 
,2 

Re* = Re, 


^ . This is the identical "reduced” or "effective” Reynolds 

number that is used in the bearing theory to show the true relative magni- 
tude of the inertia and viscous forces; therefore, 

Re* = I ner *i a forces 
Viscous forces 

(e. g. , see ref. 32). Now for curved geometries, a Reynolds number 
based on clearance has appeared. This is especially true in viscoseal 
analysis where 


Re c = 


Uc 
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Writing the previous x-direction momentum equation in terms of Re ( 
gives 


3u_ 

1 

— +— 1 
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( L 

3 2 u + c 

1 

s 2 ^ 

dz 
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0X Uc’ 
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u3h + v^; + w!h=-. 


!E + J Le. 


9 2 u , L c 
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0x 0y 0z „_2 m 0x Re„ \L 2 ~ 2 c a -2 L .,2 a -2 

cos ck c \ cos a? ox 0y sin a dz 

O _o 

Again it is seen that the dni/9y term is the largest. 

Thus l/Re c (L/c) should be of unit order. Now 


Re* = Re. 


(r) 


which leads to 


Re* = Re, 




It does not matter which Reynolds number is used to find the relative 
inertia force to viscous force ratio provided that c/L < 1, which corre- 
sponds to the conditions 


» 


and 


ay 

,2 


0z‘ 
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ay ax* 

The proper x-direction momentum equation then becomes 

0U 0U 0U 1 0P 0^U 

u + V — — + W — = — + p - 

0x 8y 0z p 0x 9y 2 


(A7) 
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It is assumed that a does not go to 0° or 90° for the formal ordering 


procedure to be valid. 


Example 


A numerical example will be used to illustrate the validity of neglect- 
ing the 3 /3z and 3 /0x viscous terms. A typical viscoseal has the 
following construction and operating conditions: smooth 2 -inch -diameter 
rotor, 0. 0035-inch concentric radial clearance between the ridge (housing) 
and rotor surfaces, 1000 rpm, and sealed fluid (water) at 60° F. There- 
fore, 

v - 17. 5X10 -4 sq in. /sec 
c = 0. 0035 in. 

U = 27tR( 1000 rpm)/60 = 100 in. /sec (plate speed) 

c/R = 0.0035 (therefore, curvature effects should be slight) 

L = Ridge width + Groove width (in circumferential direction) 

= 0. 2 in. + 0. 2 in. = 0.4 in. 
a = 15° 

Re L = UL /v = 23, 000 
Re c = Uc/ v = 200 

Using the Reynolds number based on the characteristic length L yields 

i2— / T \ 2 ^ 2— , -.2^ 


Re, 


1 3fu + /L\* cru + 1 8 u 

cos 2 o? 3x 2 c 3y^ sin^a 3z^ 


= -1- (l. 1 + 13 000 — + 15 — ^ 

Re L V 3x 2 3y 2 3z 2 / 


Using the Reynolds number based on the clearance c gives 
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Re. 


ic\ i dh + /l\ dh ] / c \ i dh 
[_ L cos^ar 9x 2 c 9y 2 ' L ' sin 2 a 9z 2 _ 

Re c V 9x 2 


2rr i2— j 

+ 114-LH+ 0. 13 £ 
9y 2 9z 


8 


0 A 

In both equations it is seen that 3 u/3y z is the largest term. Now the 
modified or reduced Reynolds number is 

.2 


Re * = Re h(i) =Re c(j;) = 1 ' 75 


Nondimensionalization of y-Direction Momentum Equation 
Nondimensionalizing the y-direction momentum equation (eq. (2)) 
yields 

u^ + v^ + wil=- f L f 3P , 1 f 1 8% , /IA 2 3^ , 1 a 2 ? 

3x By 3z Vc ) 9y Re L [ cos 2 a , g 2 U ) a _2 si „2 a 3 -2_ 

(A8) 

Again using the definition of Re* and examining the terms of unit order 


uii, viz, w — , JL 12 

9x 9y 9z Re* g -2 


a-d) 


thus 


/L\ 2 9P a 
Vc / 9y 


tf(l) 


Since 


9P 

3y 


O' 


te) 2 ] 


therefore 



f 
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— a 0 (A9) 

9y 

This result is expected for a flat plate where there is no mechanism to 
generate a radial pressure gradient (e. g. , a centrifugal force). For 
small curvatures, c/R < 0. 1 however, the radial pressure gradient 
3P/3y should also be a small quantity and have little physical inf luence 
on the axial pressure gradient. This can be shown by considering the 
following simplified model. Consider the flat plate model accelerating 
steadily about an origin. From elementary dynamics 

a = 4 t S t + a n « n 

where: a^ = q = 0 for steady acceleration 

2 2 
„ _ q _ (u cos a - w sin a) 

n R R 

q is the fluid velocity in the tangential plane (x-z plane) and 
q = u cos a - w sin a (see fig. 17). § n is the unit normal that is directed 
towards the axis of rotation. This is the y-coordinate direction in the flat 
plate model. Thus the centrifugal inertia force effects the y-momentum 
equation only. 

Add pa n to the inertia force in equation (A2) and nondimensionalize 
as before and let U= UU and R = RR. Thus equation (A8) becomes 


u — + v — + w — +-L/— ^ (- sin aw + cos au) 2 

3x 3y 3z rVc/ \R/ 



3P 

3y 


Re, 


i ^2- 

1 3 v 

2 —2 
cos a 3x 


2 -. 2 - 


/L\ 3 v 
3y2 


2_ 

1 3 v 

sin 2 a 3z 2 


(A8a) 
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Now 


2 

(— j (-sin of + cos au) 2 ~ ^(l) 


Comparing the largest viscous term, pressure and centrifugal inertia term 

^ 2 - 


ap l a 4 *? l /c \ , - _,2 

— = + — (— ] (-sin aw + cos au) 

3v Re T 3-2 r \R/ 


Now use a typical Re L (as in example) = 10 3 and apply restriction 

- < 0.01 

R~ 

~ = (10“ 3 ) — + 112- 2 ) (-sin aw + cos au) 2 
9 y ay 2 R 


Again 


3P 

3y 


“ o 


if 


£ < 0.01 

R~ 

which substantiates the claim that the centrifugal effects should be small. 
In this model Coriolis forces were also thought to be negligible. The 
above centrifugal effect is for the case when the grooves are on the sta- 
tionary housing. A slight modification is required for the groove -ridge 
pattern rotating. Note that as c/R — 1 the centrifugal force effects are 
no longer negligible. 

Nondimensionalization of z -Direction Momentum Equation 
hi a like manner the z -direction momentum equation is nondimension- 
alized and appears in a similar form as the nondimensional x-direction 
momentum equation (A6) 
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-3w + -0w + -0w = 


as 


3y 


0Z 


1 3P f 

V* 

1 

3 V 

/L\ 2 3 2 w + 1 

B 2 ^ 

. 2 3z 
sm a 

UL 

cos 2 « 

& 

' By 2 sin 2 ar 

n 2 ] 


Examining the viscous terms 

i aV 


(A 10) 


Rp. 

ii 


2 2 

cos a ox 


f L\ 2 3% + X bV I 

' c ' By 2 sin 2 u 3Z 2 ! 


The largest term is (L/c) 2 (B 2 ^/ 3y 2 ) and it would appear that the d^w/dz 2 
term could be neglected; however, this term must be retained due to a 
physical argument. Because of the "edge” there is a sharp gradient in the 
z-direction and the nature of the flow is such that the S^/Bz 2 term must 
be retained. Thus the z-direction momentum equation in dimensional form is 


3w 3w 3w 
u — + v — + w — — 

3x 


= -! 3P + „/A + 3M 

p 02 \ 3v 2 3z 2 / 


(All) 


3y 3z p az 3z < 

This was pointed out by Dr. A. Mager of Aerospace Corporation, Los 
Angeles, California (personal communication, October 5, 1966). 

Nondimensionalization of the Energy Equation 
The energy equation for steady temperature variation for an incom- 
pressible fluid with constant thermophysical properties and negligible 
compression work is 

i2m -»2r 


3T 3T 

U + V — + w 

dX 


3z pgC V ^2 ~ 2 ~2 


(A 12) 


3y 9z P% C \ax 2 3y 2 3Z 2 ) P& PZ Q 
Nondimensionalize as before and nondimensionalize the temperature by 


letting 


T - T 


e = 


s 


T» yr - T_ 
M S 


and <&' = — 
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Thus 


-a e . - ae - de 


u 


3x 


+ v 



(jL_ 

\ 

/C \ 2 l a 2 e + a 2 e + /c 

f 1 

2l 

a z e 

az 

\pgc 

Uc 2 / 

cos 2 a 9x 2 3y 2 '*■ 

‘ sin 2 a 

az 2 . 


& 


Q 


Lpg c ( T M ~ Te) pgC(T^-T c ) 


M 


now 


,PgC v 


v L\_ 1 1 /L\_ 1 = J_ 

Uc c / Pr Re„ Vc / PrRe* P* 
/ c e 


where P* is the modified Peclet number 
e 


P* = P c 
e e 


HY Ptr <z) 


or 


- ae - ae — ae 

U + V — + w — 


3X 


3y 


3z P 


(- 

f 1 

a 2 e a 2 e /c v 

l 

a 2 el 

\ L 

' cos 2 a 

* * k) 

. 2 
sm a 

az 2 _ 


$ 


Q" 


LpgC(T M -T s ) PgC(T M -T s ) 

Now assume constant temperature in the x-direction (along the groove 
ridge pairs) then 

r rt n 

0’ Q* 


v — + w 
ay 


ae _ 

i 

a 2 e 

1 ( c 

■\ 2 a 2 e 

az 

K 


sin 2 <* Vl 

az 2 . 


tpgC(T M -T s ) pgC(T M - T s ) 
Now dimensionalize all terms but the temperature 

v i£ + wi£ = -^/^ + !^V — * + — 

3y dz pgC Vay 2 dz 2 / P« c P« c 
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(A13) 


Pr =Q± 
k 

therefore 

vM +w H=j^/i!» + i!»V-K-*+-aI 

9y 0z Pr\^ 0y 2 0Z 2 J PgC pgC 

APPENDIX B 

FORMULATION OF ALONG GROOVE -RIDGE BOUNDARY CONDITION 
FOR A CYLINDRICAL GEOMETRY (VISCOSEAL) 

A special case of parallel groove geometry is a cylindrical geometry 
which is the helically grooved fluid film seal (viscoseal). For this cylin- 
drical geometry a special boundary condition must be placed on the pres- 
sure gradient along the groove-ridge direction. This can best be illus- 
trated by referring to figure 26. It is seen that the point o to point a can 
be traversed along both coordinate axes on the surface (x and z axes). 

From figure 26 and figure 21, it is seen that going from point o to 
point a along the z-axis involves the traversing of an angle cp. 


tan a = 


_ R<p 


(Bl) 


Now going the same distance along the x-axis 


Equating (Bl) and (B2) 


or 


cot a = 2lrR - 
L' 


2tt - q? _ (p 
tan a cot a 


(B2) 


<p = 2v cos a 


(B3) 
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The arc distance (path length) in the x-direction is: 

S - V (R?) 2 + < L ') 2 
Upon substitution for L* and cp 

S = 2irR cos a (B4) 

The pressure differential at a from traversing the z -direction from 
o is 

<*Va = < iP r + g)z< N ) 

then 

3T) N(AP ) N(AP ) 
oP _ v r+g'z _ v r+g 7 z 

8x S 2irR cos a 

It is seen that 9P/ dx = constant which follows from fully developed 
flow along the groove -ridge direction and verified by experiment in refer- 
ence 5. 

APPENDIX C 

SUMMARY OF CREEPING FLOW ANALYSIS 
The optimization from creeping flow analysis has been formulated by 
considering a flat -plate model (unwrapped cylinders, see fig. 27) and the 
following: 

(1) General incompressible plane Couette flow is assumed in the 
£ -direction (Poiseuille flow and simple Couette flow). The | -direction 
momentum equation is 

_ 1 3P 

gj»2 (1 0| 

with the boundary conditions 


! 



where 


h = c at ridge surface 
h = c + h 0 at groove root surface 
The well known solution (e. g. , see ref. 32) is 

-.Su-JLSEtM) 

e h 2/i d( \ hi 

(2) Poiseuille flow is assumed in the ^-direction. The 77 -direction 
momentum equation is 


0^2 71 977 

with the boundary conditions 

w =0 at £ = 0 

7 ] 

w = 0 at £ = h 
77 

The well-known solution (see ref. 32) is 

W T, = — (— 

^ 2/i. V077 / 

Now the condition Q = 0 is imposed in the control volume; that is, 

Q = + % = ff \ *> d ? + // d ? d ? = 0 

which means that there is no net end leakage or closed channel flow. This 

is a restriction that limits the optimization and pressure gradient predic- 
tion to a seal only. 

Solving for pressure gradient that satisfies the Reynolds equation 
(ref. 34) results in a sealing parameter 
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AP _ 6<nUG 

L’ 2 
c 

or 

G = -«i (**\ = I 

6jLtUVL T / A 

where A is the bearing number (ref. 34) arid G is a geometry factor. 

See references 1 to 26 for various expressions and values for G. 

The creeping flow analysis has some paradoxes. It can best be illus- 
trated by referring to figures 17 and 28 and by observing the resolution of 
the drag force into components along the groove (x -direction) and normal 
to the groove (z -direction), hi creeping flow, the drag force is in equilib- 
rium with the pressure force. Thus along the groove (x-direction) the 
drag force is indeed in equilibrium with the pressure force. In the normal 
groove direction (z -direction), however, the pressure and drag forces are 
in the same direction, which is also the direction of a very large pressure 
gradient (see figs. 22 and 23). Therefore, an unbalance of forces exists 
in this direction. Physically, the Q = 0 restriction is not possible. Also, 
this point can be illustrated from |- and 17 -direction momentum equations, 
respectively: 

d 2u e _ 1_ 9P 

0^2 

3 S _ 1 0P 

0^2 fi drq 

Notice that there is no velocity "coupling” between the velocity compo- 
nents; that is, u £ and w^ are independent of one another. It is obvious 
that for zero leakage these velocity components are not independent of one 
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another. Thus, physically as well as mathematically it is not possible 
under the restrictions of the classical creeping flow analysis to get zero 
leakage. 

Also in the creeping flow analysis the edge effects are assumed neg- 
ligible. The full meaning of this assumption has not really been clarified; 
it is only close to physical reality when c is very small and/or b/h Q is 
very large. For example, the creeping flow analysis is valid for spiral 
groove thrust bearings where c is on the order of 500 microinches and 
the aspect ratio is on the order of 1000 (c/L « 1, see fig. 29). As 
b — h Q , the edge effect becomes more pronounced as the convective forces 
becomes more and more important. To date only reference 35 has ex- 
plicitly stated this restriction. 

To resolve the paradox of having the drag and pressure forces un- 
balanced in the z -direction, a convective inertia force is in equilibrium 
with the pressure and drag forces (a nonlinear effect). The step is causing 
a convective change that results in generation of an axial pressure gradient. 
The role of the step is now described. 

If a concentric rotor is rotating about a stationary smooth sleeve, no 
axial pressure gradient is generated. This is the classic axisymmetric 
Couette flow (see ref. 36). Furthermore, if an axial pressure gradient is 
imposed (e. g. , by having the rotor translate in the axial direction) an 
’’uncoupled" helical flow will result. Subsequently, of course, there will 
be net flow out. 

In the helical groove seal (viscoseal), the step or edge h Q is the 
pressure generation mechanism. Without the step, that is, the limiting 
case when h Q — o (rotor and housing smooth), there is no axial pressure 
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gradient generated. Eccentricity will generate a radial pressure gradient. 
Therefore, the effect of the groove wall is a convective accleration or in- 
ertia effect (nonlinear effect). To neglect the convective force terms 
means that the primary axial pressure generating mechanism is neglected. 

The creeping flow solution equations are linear and pressure is a po- 
tential function. The creeping motion or slow viscous flow momentum 
equation in vector notation is 

VP = iiV 2 V 

The incompressible continuity equation says 

V • V = 0 

Taking the divergence of both sides of the equation yields 

V • VP = /iV 2 (V • V) = 0 
or 

V 2 P = 0 

which is Laplace’s equation, since pressure is a scalar quantity. Conse- 
quently, the multitude of mathematical analog methods of solving linear 
equations and specifically Laplace’s equation can be used. Thus, solu- 
tions for end effect can readily be found (refs. 1 and 34). 

Muijderman (ref. 1) used conformal mapping to analyze the spiral 
grooved bearing. Again, this technique can be readily applied because 
pressure is a potential function. 

REFERENCES 

1. Muijderman, E. A.: Spiral Groove Bearings. Philips Res. Rept. 
Suppl. 1964, no. 2. : See also Scientific Lubr. , vol. 17, no. 1, 

Jan. 1965, pp. 12-17. 



54 


2. Rowell, H. S.; and Finlayson, D. : Screw Viscosity Pumps. Engi- 

neering, vol. 114, Nov. 1922, pp. 606-607. 

3. Boon, E. F.; and Tal, S. E.: Hydrodynamic Seal for Rotating Shafts. 

DEG. Inf. Ser. 13, United Kingdom Atomic Energy Authority, 1961. 

4. King, Alan E. : Screw Type Shaft Seals for Potassium Lubricated 

Generators. IEEE Trans, on Aerospace, vol. AS-3, Supplement, 
June 1965, pp. 471-479. 

5. Ludwig, Lawrence P.; Strom, Thomas N. ; and Allen, Gordon P. : 

Gas Ingestion and Sealing Capacity of Helical Groove Fluid Film 
Seal (Viscoseal) Using Sodium and Water as Sealed Fluids. NASA 
TND-3348, 1966. 

6. Ludwig, Lawrence P.; Strom, Thomas N.; and Allen, Gordon P.: 

Experimental Study of End Effect and Pressure Patterns in Helical 
Groove Fluid Film Seal (Viscoseal). NASA TN D-3096, 1965, 

7. Vohr, J. H. ; and Chow, C. Y. : Characteristics of Herringbone- 

Grooved Gas -Lubricated Journal Bearings. J. Basic Eng. , vol. 87, 
no. 3, Sept. 1965, pp. 568-578. 

8. Billet, A. B. : Hydraulic Sealing in Space Environments. Proceedings 

of the Second International Conference on Fluid Sealing, B. S. Nau, 

H. S. Stephens, andD. E. Turnbull, eds. , British Hydromechanics 
Research Association, Harlow, Essex, England, 1964, pp. C2-17 - 
C2-36. 

9. Rowell, H. S. ; and Finlayson, D. : Screw Viscosity Pumps. Engi- 

neering, vol. 126, Aug. 1928, pp. 249-250. 



55 


10. Rowell, H. S. ; and Finlayson, D. : Screw Viscosity Pumps. Engi- 

neering, vol. 126, Sept. 1928, pp. 385-387. 

11. Rogowsky, Z. : Mechanical Principles of the Screw Extrusion Ma- 

chine. Engineering, vol. 162, no. 4213, Oct. 11, 1946, pp. 358- 
360. 

12. Strub, R. A. : Spindle Drag Pump. Machine Design, vol. 25, July 

1953, pp. 149-151. 

13. Pigott, W. T. : Pressures Developed by Viscous Materials in the 

Screw Extrusion Machine. ASME Trans. , vol. 73, Oct. 1951, 
pp. 947-955. 

14. Anon.: Theory of Extrusion. Ind. Eng. Chem., vol. 45, no. 5, 

May 1953, pp. 969-993. 

15. Eccher, Silvio; and Valentinotti, Aldo: Experimental Determination 

of Velocity Profiles in an Extruder Screw. Ind. Eng. Chem. , vol. 
50, no. 5, May 1958, pp. 829-836. 

16. Griffith, R. M. : Fully Developed Flow in Screw Extruders. Ind. 

Eng. Chem. Fundamentals, vol. 1, no. 3, Aug. 1962, pp. 180-187. 

17. Squires, P. H. : Screw Extrusion - Flow Patterns and Recent Theo- 

retical Developments. SPE Trans. , vol. 4, no. 1, Jan. 1964, pp. 
7-16. 

18. McGrew, J. M. ; and McHugh, J. D. : Analysis and Test of the Screw 

Seal in Laminar and Turbulent Operation. J. Basic Eng. , vol. 87, 
no. 1, Mar. 1965, pp. 153-162. 

19. Zotov, V. A. : Research on Helical Groove Seals. Russ. Eng. J. , 

vol. 10, Oct. 1959, pp. 3-7. 



«• 


56 

20. Asanuma, T. : Studies on the Sealing Action of Viscous Fluids. 

Paper No. A3 presented at the First International Conference on 
Fluid Sealing. Cranfield, England. British Hydromechanics Re- 
search Association, Harlow, Essex, England, Apr. 1961. 

21. Lessley, R. L. ; and Hodgson, J. N. : Low-Leakage Dynamic Seal- 

to-Space. Paper 65-GTP-14, ASME, Feb. 1965. 

22. Whipple, R. T. P. : The Inclined Groove Bearing. Rep. AERE-T/ 

R-622 (rev. ), Research Group, Atomic Energy Research Establish- 
ment, United Kingdom Atomic Energy Authority, 1958. 

23. Gruntz, Robert D. ; and Rackley, Ray A. : Snap 50/Spur Power Con- 

version System -Objectives, Current Status and Lunar Applications. 
Paper No. 650321, SAE, May 1965. 

24. Golubiev, A. I. : Studies on Seals for Rotating Shafts of High-Pressure 

Pumps. Wear, vol. 8, no. 4, July/Aug. 1965, pp. 270-288. 

25. Ho lan, Karel: Sealing in Engineering. Proceedings of the Second 

International Conference on Fluid Sealings, B. S. Nau, H. S. 

Stephens, and D. E. Turnbull, eds. , British Hydromechanics Re- 
search Association, Harlow, Essex, England, 1964, pp. E5-73 - 
E5-88. 

26. Stair, William K. : Analysis of the Visco Seal. Rep. No. ME 65-587-2, 

University of Tennessee, Jan. 18, 1965. 

27. Stair, William K. : Theoretical and Experimental Studies of Visco- 

Type Shaft Seals. Rep. No. ME 66-587-5, University of Tennessee, 
Apr. 28, 1966. 



57 


28. Lessley, R. L. : Snap-8 Seals-to -Space Development Test Program. 

Vol. 1 - Visco Pump. Rep. No. 2808 (Topical) (NASA CR-54234), 
Aerojet-General Corp. , May 1964. 

29. Kettleborough, C. F.: Turbulent and Inertia Flow in Slider Bearings, 

ASLE Trans. , vol. 8, no. 3, July 1965, pp. 286-295. 

30. Zuk, John; Ludwig, Lawrence P. ; and Johnson, Robert L. : Flow and 

Pressure Field Analysis of Parallel Groove Geometry for an Incom- 
pressible Fluid With Convective Inertia Effects. NASA TN D-3635, 
1966. 

31. Kuivinen, David E.: Determination of Oxygen in Liquid Alkali Metals 

by the Mercury Amalgamation Method. Paper presented at the 
Nineteenth Meeting, Chemical Rocket Propulsion Group (St. Paul, 
Minn. ), July -Aug. , 1963. 

32. Schlichting, Hermann (J. Kestin, Trans. ): Boundary Layer Theory. 

Fourth ed. , McGraw-Hill Book Co. , Inc. , 1960. 

33. Kulinski, E. S. ; and Ostrach, S.: A Critical Evaluation of the Status 

and Trends in High Speed Fluid Film Lubrication (to be published). 

34. Bisson, Edmond E. ; and Anderson, William J. : Advanced Bearing 

Technology. NASA SP-38, 1964. 

35. Booy, M. L. : Influence of Oblique Channel Ends on Screw -Pump 

Performance. J. Basic Eng. , vol. 88, no. 1, Mar. 1966, pp. 
121-131. 

36. Langlois, W. E. : Slow Viscous Flow. Macmillan Co. , 1964, pp. 


105-107. 



V 


TABLE I. - VISCOSITY 


OF MINERAL OIL 
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^Used as base line. 
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TABLE IH. - CRITICAL MODIFIED REYNOLDS NUMBER FOR ONSET OF TURBULENCE OF VARIOUS EXPERIMENTERS 


Experimenter 

Reference 
and fluid 

Seal 

Diameter, 

in. 

Radial 

clearance, 

c, 

in. 

Ratio of 
radial 

clearance to 
rotor radius, 
c/R 

Ridge 

width, 

a, 

in. 

— 

Groove 

width, 

b, 

in. 

Helix 

angle, 

a, 

deg 

Characteristic 

length, 

L, 

in. 

— 

Clearance 

Reynolds 

number, 

Modified 

Reynolds 

number. 

Re* 

Stair 

27 

a l 

1.2430 

0.0042 

0.0068 

0.0934 

0.1596 

14.5 

1.01 

400 

1.7 


Water 

2 

1.2465 

.00235 

.0038 

.1176 

.0494 

9.67 

.98 

(b) 

(b) 



2B 

1.2420 

.0047 

.0076 

.1176 

.0494 

9.67 

.98 

(b) 

Cb) 



3 

1.2461 

.00265 

.0043 

.0828 

.0842 

9.67 

.98 

350 

.95 



3B 

1.2408 

.0053 

.0085 

.0828 

.0842 

9.67 

.98 

(b) 

<b> 



4 

1.2461 

.00265 

.0043 

.0514 

.1156 

9.67 

.98 

200 

.95 



4B 

1.2408 

.0053 

.0085 

.0514 

.1156 

9.67 

.98 

400 

2.2 



5 

1.2455 

.00295 

.0047 

.1585 

.0607 

5. 81 

1.83 

600 

.97 



6 

1.2455 

.00295 

.0047 

.1070 

.1083 

5. 81 

1.83 

200 

.32 



7 

1.2455 

.00295 

.0047 

.0609 

.1600 

5. 81 

1.83 

160 

.26 

King 

4 and 27 

3C 

2 

0.0029 

0.0029 

0. 1063 

0. 1063 

3.86 

3.2 

550 

0.50 


Oil, 

IB 

2 

.0020 

.0020 

.1063 

.1063 

3.86 

3.2 

400 

.25 


water, po- 

2E 

2 

.0025 

.0020 

.1000 

.1000 

7. 26 

1.54 

650 

1.1 

| 

tassium d 

5D 

2 

.0015 

.0015 

.1000 

.1000 

7.26 

1.54 

650 

.63 

McGrew and 

18 

SF96-5 oil 

1 

0. 00325 

0.0065 

(c) 

(c) 

14.5 

3. 14 

400 

0.41 

McHugh 

Oil 

j 










Lessley 

28 

Configura- 

2 

0.0032 

0. 0032 

0. 075 

0.122 

14.5 

0.79 

350 

1.4 


Oil 

tion B 











Water 

Configura- 

2 

.0032 

.0032 

.083 

.136 

4 

3.13 

500 

.51 


Mercurj^ 

tion C 











a Eccentricity of rotor to bousing, c = 0. 1. 

^Break point not well defined. 
c Four thread starts. 

^Not explicitly found at breakpoint but falls on data line. 



HJ-3by± 


# « 
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| Water out 



| To leak detector 
r* 

Liquid-gas interface-^ 
VA y Argon gas ' 


-Land width, a* 


Groove width, b* 


Argon gas 
atmosphere 
t 



Helix angle, a 


Line to reservoir 


Figure 4. - Grooved rotor viscoseal used to seal sodium. 
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Figure 5. - Sodium supply system. 



Pressure/Cavfty pressure 


CD 

to 

K) 

i 

W 





Circumferential distance, rad 


(b) Rotational speed, 4000 rpm; average temperature of oil film, 138° F; temperature change along 
film, 4° F; cavity pressure, 60. 2lb/sqin. absolute. 

Figure 7. - Pressure gradients on planes orthogonal to viscoseal axis centerline. 
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Pressurized __ 
cavity 

End of seal 




3 



(a) Rotational speed, 1000 rpm; average temperature of oil film, 
97° F; temperature change along oil film, 11° F; cavity pres- 
sure, 54. 1 Ib/sq in. absolute. 



Axial distance 


(b) Rotational speed, 6000 rpm; average temperature of oil film, 
132° F; temperature change along oil film, 19° F; cavity pres- 
sure, 80. 6 Ib/sq in. absolute. 

Figure 8. - Axial pressure gradients near helical groove edges 
with grooves in direct communication with pressurized 
cavity. 
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(a) Modified Reynolds number, 1.75 (2000 rpm). 
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(b) Modified Reynolds number, 2.58 ( 3000 rpm). 

Figure 9. - Grooved housing viscoseal operation when sealing water. 
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(c) Modified Reynolds number, 3.50 (4000 rpm). 
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(d) Modified Reynolds number, 4.42 (5000 rpm). 


Figure 9. - Continued. 
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(e) Modified Reynolds number, 5.25 (6000 rpm). 
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(f) Modified Reynolds number, 7.00 ( 8000 rpm). 



Figure 9. - Concluded. 
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Modified Reynolds number 

Figure 10. - Comparison of observed gas ingestion 
when sealing water for grooved housing and 
grooved rotor viscoseals. 



(a) Grooved housing. Bubble motion pri- 
marily circumferential. 

Figure 11. - Gas ingestion mechanisms of grooved 



(b) Grooved rotor. Bubble motion has axial 
component due to helical grooves. 

housing compared with grooved rotor viscoseal. 
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Figure 12. - Pressure as function of axial length wi th gas ingestion for 
grooved rotor. Sealed fluid, water. 
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Figure 13. - Gas ingestion rate when sealing liquid sodium at 400° F. 
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(see fig. 5) 

(a) Mechanism involved in gas ingestion into pressurized cavity. 



(b) Schematic liquid level oscillograph trace for repeated cycles of gas ingestion, subsequent 
gas blowout, and no liquid leakage at blowout. 



(c) Schematic liquid level oscillograph trace for repeated cycles of gas ingestion, subsequent 
gas blowout, and liquid leakage at blowout. 

Figure 14. - Gas ingestion into closed cavity. 


E-3691 



10 3 10 4 10 5 


Reynolds number based on clearance 

Figure 15. - Comparison of sealing parameters for grooved rotor and housing 
operating in sodium. Sodium temperature, 165° to 335° C (329° to 635° F). 
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(b) Cross-sectional view along z-axis illustrating relative 
motion of smooth flat plate with respect to fixed parallel 
groove-ridge plate. 

Figure 17. - Mathematical model of parallel groove- ridge geometry of infinite horizontal extent 
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Figure 18. - Region of analysis validity (homogeneous fluid) of a grooved housing viscoseal. Reynolds 
number based on clearance or film thickness, 1700; modified or reduced Reynolds number, 19; 
400-microinch shaft rotational movement during photographic exposure; shaft speed, 8000 rpm : 
sealed fluid, water. 
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Figure 19. - Viscoseal with grooved housing operating in a regime of water-air mix- 
ture (gas ingestion). Analyses does not apply for this discontinuous fluid film 
case; 250-microinch shaft rotational movement during photographic exposure; 
shaft speed, 5000 rpm; sealed fluid, water. 



Figure 20. - Cross-sectional view along z-axis 
illustrating arbitrary groove-ridge pair shape. 


Circumferential 

Figure 21. - Resolution of across and along groove-ridge pressure gradients into 
axial pressure gradient "Unwrapped 11 cylinder is shown. 



Figure 22. - Pressure distribution over parallel groove geometry. (Qualitative 
graphical representation from experimental data.) 
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(b) Across ridges and grooves (no restric- 
tion that profile must be linear). 

Figure 23. - Pressure profiles (qualitative graphical representation from experimental data). 
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Clearance Reynolds number 

Figure 24. - Experimental sealing parameter for seals 5, 6, and 7 from reference 27. 
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Figure 29. - Example where the creeping flow mathematical model can be considered a valid 
physical model. Modified Reynolds number (Re*) much less than 1 since clearance c is 
much less than characteristic length L. 
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